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MINIMUM RADII OF SUPER-EARTHS: CONSTRAINTS FROM GIANT IMPACTS 
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ABSTRACT 

The detailed interior structure models of super-Earth planets show that there is degeneracy in 
the possible bulk compositions of a super-Earth at a given mass and radius, determined via radial 
velocity and transit measurements, respectively. In addition, the upper and lower envelopes in the 
mass-radius relationship, corresponding to pure ice planets and pure iron planets, respectively, are 
not astrophysically well motivated with regard to the physical processes involved in planet formation. 
Here we apply the results of numerical simulations of giant impacts to constrain the lower bound in 
the mass-radius diagram that could arise from collisional mantle stripping of differentiated rocky/iron 
planets. We provide a very conservative estimate for the minimum radius boundary for the entire 
mass range of large terrestrial planets. This envelope is a readily testable prediction for the population 
of planets to be discovered by the Kepler mission. 

Subject headings: planets and satellites: formation — planetary systems 



1. INTRODUCTION 

Since the first confirmed detection of an extrasolar 
planet in 1995, the catalog of known exoplanets has 
grown to surpass 400. Precision Doppler shift discovery 
techniques have now identified more than a dozen exo- 
planets in the mass range from about 2 to 10-15 Earth 
masses (M^). They are expected to be terrestrial in na- 
ture and unlike the gas giant planets in our solar system. 
This new class of planets has been collectively termed 
super- Earths" (jMelnick et al.ll2001[ ). Recently two tran- 
siting super-Earths were discovered, providing for the 
first time r adii, in addition t o masses. One of th em — 
CoRoT-7b (jLeger et all [200l iQueloz et all [2009( )— has 
high density and is likely ro cky. The other one — GJ1214b 
(jCharbonneau et al.ll2009l ) — has low density and is likely 
water rich and surrounded by a small hydrogen-rich en- 
velope. Many more super-Earths and measurements 
of their radii are e xpected from the Kepler mission 
(jBorucki et al.l[200l . 

Theorists anticipate a rich diversity in the bulk com- 
position and internal structure of super-Earths, which 
is reflected in the broad band of possible radii on a 
planetary mass-radius diagram. The band corresponds 
to the anticipated range of mean densities. There are 
four distinct types of materials that could make up a 
planet in this regime: silicates, iron alloys, volatiles/ices, 
and hydrogen-helium gas. The range of possible mixing 
ratios between these materials leads to degeneracies in 
the determi nation of bulk composition from radius and 
mass alon e (ISasselov et al 1120081: lAdams et al.ll2008D . As 
shown by [Valencia et al.l (|2007[ ). in order to restrict the 
range of possible bulk compositions, precise radii and 
masses (to 5% and 10%, respectively) have to be comple- 
mented with knowledge of stellar abundance ratios (e.g., 
Si/Fe) and physical constraints on the maximum fraction 
of H2O or iron in a planet. The latter constraints place 
limits on the maximum and minimum possible radii for 
solid planets, respectively. 
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In this Letter, we consider the minimum possible ra- 
dius a super-Earth could have at a given mass. Since 
we are interested in the limiting case, our discussion can 
be confined to rocky planets composed of iron and sil- 
icates with no ices/water or hydrogen-helium gas lay- 
ers. Under physically plausible conditions around normal 
stars, planet formation will lead to differentiated super- 
Earths with an iron core and a silicate mantle, with the 
proportions of each determined by the local Si/Fe ratio 
(jCrasset et al.l [20091 ). The only way to significantly in- 
crease the mean density requires removal of the silicate 
mantle while preserving the iron core. The most efficient 
method to strip the mantle is by giant impacts, which 
are common in the final stages of planet formation (e.g., 
IChamberil2004f) . Given the large gravitational potential 
of super-Earth planets, we suspect that complete mantle 
stripping is not possible, and unlike the case of aster- 
oids, pure iron super-Earths do not exist around normal 
stars. In this Letter, we analyze the results from numeri- 
cal simulations of planet-planet collisions and determine 
a theoretical lower limit on the planetary radii of super- 
Earths. We anticipate that observations by the Kepler 
mission will test our predictions. 

The initial conditions for our investigation are depen- 
dent on an understanding of planet formation. These 
results and observations could, in turn, help constrain 
theories fo r planet formation. In re cent years. ITda fc LinI 
(2004) and lMordasini et al.l (|2009al ) have applied detailed 
models of planet formation to the generation of syn- 
thetic populations. Such synthetic populations can then 
be compared to the observed distribution of known ex- 
oplanets to place statistically significant constraints on 
planet formation model s. For example, the analysis of 
iMordasini et al.l (|2009ai rbl) shows that the core accretion 
model of planet formation can reproduce observed popu- 
lations. In this model, small planetesimals collide to form 
larger planetary embryos O.IM0), the most massive of 
which come to dominate accretion in a process known as 
runaway growth. This stage is followed by an oligarchic 
stage in which protoplanets become relatively isolated 
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after consuming the surrounding planetesimals. At this 
stage, cohisions between comparably sized large bodies, 
giant impacts, become important and dominate the end 
stages of the formation of terrestrial planets. Giant im- 
pacts have been extensively modeled in the planetary 
embr yo size regime (jAgnor &: A sphaug 2004; lAsphau g 
120091 ): however they had not been studied extensively up 
to Earth size , with the ex ceptions of the Moon forming 
impact (e.g.. ICanud 120041 ). Recently, the first study fo- 
cused on collisions between super-Earths determined the 
criteria for catastrophic disruption and deriv ed a scaling 
law for mantle stripping (|Marcus eTaDEOOi). 

2. METHOD AND ASSUMPTIONS 

To address the question of physically plausible mini- 
mum radii for super-Earths, we consider (1) the mini- 
mum radius based on cosmic chemical abundances and 
(2) subsequent mantle stripping by giant impacts. 

The initial state of super-Earths is ba sed on two as- 
sumptions, following the discussion in iValencia et al.l 
dlOOT). The first is that all super-Earth sized planets 
have undergone differentiation. The second is that the 
relative elemental abundances are the same as those in 
the solar system and the solar neighborhood. The first 
assumption is not restrictive since all the terrestrial plan- 
ets and large satellites in the solar system are known 
to be differentiated. Further, if any super-Earth planet 
were not fully differentiated at the time of the impact 
considered here, it would be impossible to preferentially 
remove lighter materials, and thus our constraint on the 
minimum radius possible from collisional mantle strip- 
ping would still hold. The second assumption should 
also be adequate for the super-Earths that we expect to 
be discovered in the near future (close-by and in the pre- 
s elected targets o f the Kepler mission) . 

IValencia et al.l (|2007| discuss the process by which 
cosmic elemental abundances constrain an initial mini- 
mum radius for super-Earths, which we summarize here. 
Volatiles, silicates and metals condense at different tem- 
peratures. As the stellar nebula cools, the most re- 
fractory elements condense first. First, silicates con- 
dense at temperatures between 1750 and 1060 K, fol- 
lowed b y the metals (e.g., Fe , Ni) between 1450 and 
1050 K (jPetaev &: WoodI 120051 ). depending on the pres- 
sure in the nebula. Finally, H2O and other ices condense. 
We estimate the maximum iron core mass fraction by 
considering the relative mass fraction of major elements 
in the solar nebula: H at 74%, O at 1.07%, Fe at 0.1%, 
Si at 0.065% and Mg at 0.058%. We consider Si and 
Mg to be practically equally abundant. During the con- 
densation sequence (for pressures < 10"** bars). Si wiU 
condense before Fe. If Fe remains immiscible, the largest 
core is attained at a mass ratio of Si/Fe~ 0.6. In this 
case, the mantle is effectively MgO+Si02 so that the 
Si/Fe ratio can be used as a proxy for the mantle-to-core 
mass fraction. 

We contemplate processes in the late stages of planet 
formation that influence the final state of a planet start- 
ing from embryos with normal cosmic abundances (an 
iron mass fraction of about 0.33, which is the value for 
Earth). Any process that preferentially induces the es- 
cape of light elements (e.g., solar wind, gravitational es- 
cape) will deplete the planet from volatiles, including 
H2O, and perhaps silicates. The major widespread pro- 



cess that could change the mantle-to-core mass fraction 
dictated by elemental abundances is giant impacts. The 
most effective collisions for increasing the mean density 
of a planet are likely to be near-equal mass planet-planet 
encounters, where both bodies are dry (composed of iron 
and silicates). There would be few such planets, es- 
pecially in the super-Earth mass range, in any planet- 
forming disk. Hence, such collisions would be rare oc- 
currences, rather than multiple in the history of a given 
planet. 

Further assumptions in our calculations include the 
following: (1) the mantle is stripped in a single, late 
giant impact, when the planet is nearly fully accreted; 
(2) almost all of the mass remaining in the post-impact 
planet is from the largest remnant (in other words, the 
smaller fragments are not re-accreted); (3) super-Earths 
form only in the mass range 1-15 and beyond this 
upper limit, runaw ay gas accretion causes the planet to 
become a gas giant (|Ida fc Linll2004D . If super-Earths do 
form at masses larger than 15 M®, the physics of colli- 
sional mantle stripping would not alter, and thus our re- 
sults could be extrapolated beyond this limit. The large 
relative velocities necessary to strip a significant portion 
of a planet's mantle are most likely to occur early in 
the final stage of planet formation, during which time 
there are many small planetary embryos (~ O.IM^). In 
A^-body simulations of terrestrial planet formation, this 
stage has been shown to produce impact velocities as 
high as six times the mutual escape velocity (|Agnor et al.l 
Il999f ). However, such high impact velocities are still rare, 
so most of the planetary embryos that are eventually in- 
corporated into a fully accreted super-Earth will not have 
suffered mantle-stripping collisions, thus erasing the sig- 
nature of a small number of such impacts. 

While relative impact velocities are expected to be 
around the mutual escape velocity when only super- 
Earth mass planets remain, there need only be a single 
impact event to leave a large remnant with a high iron 
mass fraction. Note that clearing the smaller fragments 
from the orbit of this largest remnant, so that they would 
not be re-accreted, would most likely require the presence 
of either the protoplanetary disk or another large planet 
in the system. Because we are interested in presenting a 
lower limit for the radii of super-Earths as a function of 
mass, we believe these assumptions are justified and will 
at worst result in an underestimation of the minimum of 
the mass-radius relationship. 

The calculations of mantle stripping presented in this 
Letter are derived from simulations of head-on collisions 
between super-Earths. Such low angle impacts are most 
effective at stripping mantle material (|Benz et al.lll988l 
2007; Marcus ct al. 2009). The impact angle has little ef- 
fect on the efficiency of mantle stripping between impact 
angles of 0° (head-on) and about 30°. In this regime, 
highly disruptive impact events are possible. As the im- 
pact angle increases beyond 30°, there is a sharp tran- 
sition, b eyond which collis i ons e n ter the "hit-and-run" 
regime (Agnor & Asphaug 2004; 'A sphaug et al.l 120061 : 
lAsp haug 2009; Marcus ct al. 2009), in which the pro- 
jectile and target emerge from the impact largely intact. 

3. RESULTS 

Using the velocity-dependent catastrophic disruption 
criteria of .Stewart fc Leinhardt (,2009.) and the scaling 
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law for changes to the iron-to-sihcate ratio via disrup- 
tion of bodies with an in itial iron mass fraction of 0.33 
froml Marcus et all ()2009D . we derive the impact velocity 
(in km s~^) necessary to produce a post-impact largest 
remnant of a specified iron mass fraction fp,,: 

T/.^10.5(M-0.33)°-^°^(^ ^ j Ml,^. (1) 

Here, 7 is the projectile-to-target mass ratio and Mtarg is 
in Earth masses. The mass of the corresponding largest 
remnant is 

Mir = [-iXfPe - 0.33)Tfe + l]Mtargil + l) ■ (2) 

These eq uations are obtained by combining Equations 
(2)-(4) of' Marcus et"aLl (120091) 

Figure [T] presents this critical velocity as a function of 
the mass of the largest remnant, which would be the ob- 
served super-Earth. In Figure [TJa), the largest remnants 
have fpe ~ 0.7, making these planets super-Mercuries. 
For projectile-to-target mass ratios < 1/4, the critical ve- 
locity rapidly exceeds 50 km s~^, making such collisions 
all but impossible in the vicinity of 1 AU (the maximum 
possible impact velocity at 1 AU around a solar-like star 
is ~ 70 km s~^). Such large impact velocities are compa- 
rable to the orbital velocity at a location of 0.1 AU (for 
a IMq star), the location at which super-Earth plan- 
ets may end Type I migration (jMasset et al.ll2006l ). For 
masses close to IM^ and projectile-to-target mass ratios 
> 1/2, the critical velocities are < 40 km s~^. In Figure 
[ijb), the projectile-to-target mass ratio is fixed at 1:1, 
the limiting (most destructive) case. The impact veloc- 
ities necessary to produce super-Earths with 40%-50% 
iron by mass are 15-35 km s~^. 

Figure [2] presents the iron mass fraction, fpe, of the 
largest remnant after an impact as a function of the 
projectile-to-mass ratio, impact velocity, and mass of the 
largest remnant. The largest possible iron mass frac- 
tion resulting from mantle stripping in a giant impact is 
clearly a function of the final mass of the planet and the 
impact conditions. As expected, it is far easier to remove 
mantle material from a lower mass super-Earth. 

Next, the minimum mass-radius relationship for super- 
Earths can be reconsidered based on the likelihood of the 
impact conditions necessary to achieve a certain bulk 
density (based on the iron mass fraction). Figure |3] 
presents the radii of the post-impact remnants as a func- 
tion of the projectile-to-mass ratio and the impact ve- 
locity. The radii were calculated from th e results of the 
super-Earth internal structure models of [Valencia et al.l 
dlOOe, 2007). Thus, these radii are not the radii of the 
largest remnants immediately after the collision, at which 
time the largest remnant consists of gravitationally reac- 
cumulated debris (a rubble pile). Rather, these radii cor- 
respond to the radius of the super-Earth with a core mass 
fraction given by fpe long after differentiation of the reac- 
cumulated body and radiative loss of the excess thermal 
energy from the collision. This corresponds to the time 
at which we are most likely to observe the planet, given 
that formation takes only ~ 10 Myr of the planet's billion 
year-plus lifetime. From Figure [31 it is clear that even 
given quite extreme impact conditions, with velocities of 
up to 80 km s~^, the minimum curve in the mass-radius 



diagram lies well above the case for pure iron (lower black 
line), particularly for more massive super-Earths. 

Figure [4] presents the mass-radius diagram for super- 
Earths, with upper and low e r env elopes (dotted lines) as 
calculated bv' Fortnev et al.l ()2007f ). We add to this a new 
lower constraint on the radius from collisional mantle 
stripping (solid line), corresponding to the lowest blue 
line in Figure[3l a head-on collision at 80 km s~^ between 
equal-mass bodies. 

4. DISCUSSION AND CONCLUSIONS 

Using the results of giant impact simulations presented 
in Marcus et al. (2009), we have described the im- 
pact conditions necessary to strip away mantle mate- 
rial from a nearly fully accreted planet. Combining this 
with detailed interior structure models fo r super-Earths 
(jValencia et al.|[2006l [20071: iFortnev et all l2007) , we have 
constructed a mass-radius diagram for super-Earths and 
shown that the previous lower envelope in this relation- 
ship, corresponding to 100% iron super-Earths, is not 
consistent with collisional mantle stripping. Further, if 
the absence of super-Earths in t he 10-lOOMm ran ge seen 
in the planet formation models of llda fc LinI ()2004[ ) is cor- 
rect, even the existence of super-Mercuries, ~ 70% iron 
by mass, may be limited to masses < 5Mq (as with the 
top blue curve in Figure [2ja)). This restriction arises be- 
cause IOMq target bodies are required to make super- 
Mercuries larger than about 5M^. 

The Kepler mission will discover a few hundred planets 
in the mass-radius range of super-Earths shown in Fig- 
ure 4. We predict that the lower envelope of the distribu- 
tion that Kepler is going to measure will be significantly 
higher than our computed minimum based on collisional 
stripping. Our calculation derives a very conservative 
limit with very low probability for these extreme sce- 
nario to be realized. Hence, Kepler's limited sample of 
planets is unlikely to be large enough to include such rare 
events. If super-Earths violating the collisional stripping 
limit are actually confirmed around normal stars, then 
we would need to revisit basic assumptions about the 
planet formation process. 

The simulations in this Letter were run on the Odyssey 
cluster supported by the Harvard FAS Research Comput- 
ing Group. 




Fic. 1. Critical impact velocity required to obtain the specified 
post-impact iron mass fraction vs. mass of the largest impact rem- 
nant, (a) Largest remnant with iron mass fraction of 70% (sim- 
ilar to Mercury) for various projectile-to-target mass ratios, (b) 
Largest remnant with various iron mass fra<;tions for equal mass 
collisions at the lowest velocity for such an enrichment. Note that 
the solid line in both panels is for a 1:1 mass ratio collision yielding 
a largest remnant that is 70% iron. 
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Fig. 2. — Iron fraction of the largest remnant vs. mass of the 
largest remnant, (a) Impact velocities given in terms of the mutual 
escape velocity, (b) Impact velocities in km s~^. In both panels, 
line color indicates the projectile-to-target mass ratio: green, 1/10; 
red, 1/4; blue, 1/1. 
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Fig. 3. — Mass-radius diagram for super-Earths, (a) Impact ve- 
locities given in terms of the mutual escape velocity, (b) Impact 
velocities given in km s~^. In both panels, the two solid black 
lines represent terrestrial composition (iron mass fraction of 0.33, 
upper line) and pure iron (lower line) (Valencia ct al. 2007). The 
colored lines are for super-Earths stripped of mantle material in 
a giant impact. Line color indicates the projectile-to-target mass 
ratio: green, 1/10; red, 1/4; blue, 1/1. In (b), the X symbols repre- 
sent potentially observable transiting planets with 5% uncertainty 
in the radii. Note that the difference in minimum radii between 
a pure iron planet and plausible densities achieved via coUisional 
stripping is observable. 
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Fig. 4. — Mass-radius diagram for super-Earths. The dotted lines 
are pure iee (upper) and pure iron (lower) (Portney et al. 200'3). 
The solid line is a eonstraint placed on the mass-radius diagram 
from coUisional mantle stripping, corresponding to a 80 km s"-*^ im- 
pact with a 1:1 projectile-to-target mass ratio. The extrapolation 
(dashed) corresponds to cases that require target masses > 15M(g. 
The X symbols indicate the masses and radii of Earth, Uranus, 
and Neptune. 



8 



REFERENCES 



Adams, E. R., Seager, S., & Elkins-Tanton, L. 2008, ApJ, 673, 
1160 

Agnor, C, & Asphaug, E. 2004, ApJ, 613, L157 
Agnor, C. B., Canup, R. M., & Levison, H. F. 1999, Icarus, 142, 
219 

Asphaug, E. 2009, Annual Review of Earth and Planetary 

Sciences, 37, 413 
Asphaug, E., Agnor, C. B., & Williams, Q. 2006, Nature, 439, 155 
Benz, W., Anic, A., Horner, J., & Whitby, J. A. 2007, Space 

Science Reviews, 132, 189 
Benz, W., Slattery, W. L., & Cameron, A. G. W. 1988, Icarus, 

74, 516 

Borucki, W.. ct al. 2009, in lAU Symposium, Vol. 253, lAU 

Symposium, 289-299 
Canup, R. M. 2004, Icarus, 168, 433 

Chambers, J. E. 2004, Earth and Planetary Science Letters, 223, 
241 

Charbonneau, D., et al. 2009, Nature, 462, 891 
Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 
1661 

Grasset, O., Schneider, J., & Sotin, C. 2009, ApJ, 693, 722 
Ida, S., & Lin, D. N. C. 2004, ApJ, 604, 388 
Leger, A., et al. 2009, A&A, 506, 287 

Marcus, R. A., Stewart, S. T., Sasselov, D., & Hernquist, L. 2009, 
ApJ, 700, L118 



Masset, F. S., Morbidelli, A., Crida, A., & Ferreira, J. 2006, ApJ, 
642, 478 

Melnick, G. J., et al. 2001, in Bulletin of the American 
Astronomical Society, Vol. 34, Bulletin of the American 
Astronomical Society, 559 
Mordasini, C, Alibert, Y., & Benz, W. 2009a, A&A, 501, 1139 
Mordasini, C, Alibert, Y., Benz, W., & Naef, D. 2009b, A&A, 
501, 1161 

Petacv, M. I., & Wood, J. A. 2005, in Astronomical Society of the 
Pacific Conference Series, Vol. 341, Chondrites and the 
Protoplanctary Disk, ed. A. N. Krot, E. R. D. Scott, & 
B. Rcipurth, 373-406 

Queloz, D., ct al. 2009, A&A, 506, 303 

Sasselov, D. D., Valencia, D., & O'ConncU, R. J. 2008, Physica 

Scripta, T130, 014035 
Stewart, S. T., & Leinhardt, Z. M. 2009, ApJ, 691, L133 
Valencia, D., O'Connell, R. J., & Sasselov, D. 2006, Icarus, 181, 

545 

Valencia, D., Sasselov, D. D., & O'Connell, R. J. 2007, ApJ, 665, 
1413 



